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We formulate a microscopic theory of non-local electron transport in three-terminal diffusive 
normal-superconducting-normal (NSN) structures with arbitrary interface transmissions. At low 
energies e we predict strong enhancement of non-local spectral conductance 1712 oc 1/e due to 
quantum interference of electrons in disordered N-terminals. In contrast, non-local resistance R12 
remains smooth at small e and, furthermore, is found to depend neither on parameters of NS 
interfaces nor on those of N-terminals. At higher temperatures R12 exhibits a peak caused by the 
trade-off between charge imbalance and Andreev reflection. Our results are in a good agreement with 
recent experimental observations and can be used for quantitative analysis of future experiments. 



In hybrid NS structures quasiparticle current flowing 
in a normal metal is inevitably converted into that of 
Cooper pairs inside a superconductor. For quasiparti- 
cle energies above the superconducting gap e > A this 
conversion is accompanied by electron-hole (or charge) 
imbalance-'^ which relaxes inside a superconductor at a 
typical inelastic length usually denoted as Ag* . As a re- 
sult, at temperatures near the critical one Tc an electric 
field penetrates into a superconductor causing resistance 
enhancement for NS structures under consideration. 

At subgap energies e < A the physical picture be- 
comes entirely difi^erent. In this case quasiparticle-to- 
Cooper-pair current conversion is provided by the mech- 
anism of Andreev reflection^: A quasiparticle enters the 
superconductor from the normal metal at a length of or- 
der of the superconducting coherence length ^5, forms a 
Cooper pair together with another quasiparticle, while a 
hole goes back into the normal metal. Due to this process 
subgap conductance of the NS structure remains non-zero 
down to T = 0^. Furthermore, in the presence of disor- 
der this subgap conductance can be greatly enhanced at 
low energies due to quantum interference effects^'^'^. 

Further interesting effects may occur in three-terminal 
NSN structures. Provided the distance between two N- 
terminals is smaller than or comparable with ^5, elec- 
trons penetrating into the superconductor from the first 
N-terminal may form Cooper pairs with electrons from 
the second N-terminal. Then a hole goes into the second 
N-metal making the charge transfer effectively non-local. 
This important phenomenon of non-local (or crossed) 
Andreev refiection (CAR)^ enables direct experimental 
realization of entanglement between electrons from spa- 
tially separated N-terminals. 

CAR was detected and investigated in several recent 
experiments^'^'^iii by measuring the non-local resistance 
of multiterminal NSN systems. The authors observed a 
rich structure of different features many of which are still 
waiting for their theoretical interpretation. Note that not 
only CAR but also other physical processes contribute to 
the non-local conductance gi2 thus making this interpre- 
tation rather complicated. For instance, the contribution 
of clastic cotunneling (EC) to gi2 exactly cancels that of 
CAR in the lowest order in NS interface transmissions 
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FIG. 1: NSN structure under consideration. 



and at subgap energiesiS.. This cancellation is lifted ei- 
ther in higher orders in barrier transmissionslS. or in the 
presence of interactions, e.g., with an effective external 
environmenlji^, or under external ac bias^^. 

Another important issue is the effect of disorder in 
metallic terminals which needs to be analyzed for ade- 
quate interpretation of experimental results^i^ '^^i^^ . Al- 
though CAR in disordered NSN structures was already 
addressed in a number of theoretical work o^^i^^'^^i^^'^° in 
various physical limits, we believe that general analysis of 
this issue is still missing in the literature. For instance, 
the role of disorder-induced electron interference^*^ in 
the non-local subgap transport, the effect of high NS 
barrier transmissions on CAR as well as some other fea- 
tures remain unclear. Yet another important unresolved 
problem is to describe an interplay between CAR and 
non-local charge imbalance. It was demonstrated both 
experimentally®'J^° and theoretically^^'^^ that such inter- 
play may result in a large non-local resistance peak which 
occurs at temperatures slightly below the critical one. It 
was conjectured^^ that the behavior of this peak is con- 
trolled by the charge imbalance length Ag* parametri- 
cally exceeding the length scale ~ relevant for CAR. 
In this work we develop a general theory of non-local 
electron transport in diffusive NSN structures which en- 
ables one to clarify the above issues and to formulate 
predictions to be tested in future experiments. 

Model and basic equations. In what follows we will an- 
alyze a multiterminal diffusive NSN structure schemat- 
ically shown in Fig. 1. Two normal terminals iVi and 
A'2 with resistances r^fj and rjv^ and electric poten- 
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tials Vi and V2 are connected to a superconducting elec- 
trode of length L with normal state (Drude) resistance 
vl and electric potential V = via small NS barriers 
with resistances Ri and R2 which can be expressed via 
channel transmissions Ti^„ and T2^„ of these barriers as 
1/^1(2) = (g^/tt) X)n ^i(2),n- For the sake of definite- 
ness in Fig. 1 we chose specific geometry directly re- 
lated, e.g., to experimentslS where the superconductor 
was fabricated in the form of a rather thin strip. NS bar- 
riers are located at the points ri 2 = (a;i^2:07 0) and the 
corresponding segments of a superconducting strip have 
normal state resistances r^^ and r^-x^i see Fig. 1. 

Our analysis is based on the quasiclassical Usadel equa- 
tions for the Green-Keldysh matrix functions G. In the 
absence of interactions these equations read^^ 



iD\/{GVG)^[t,G\, G^^l, 
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where [a, 6] = a6 — ha, D is the diffusion constant and 
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are 4x4 matrices in Keldysh(g)Nambu space, e is the 
quasiparticle energy, A(T) is the superconducting order 
parameter which will be considered real further below 
and V is the electric potential. 

Far from the interfaces between metals the quasiclas- 
sical Green functions G coincide with their bulk equilib- 
rium values. Deep in the superconductor they read 
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Gf = (Gf -G|M£), (3) 



where 71(e) = tanh(e/2T) and fi are Pauli matrices. In 
the normal terminals far from the tunnel barriers one has 
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while the retarded and advanced Green functions G^2'^ 
are set by the first Eq. ([3]) with A = 0. In the vicin- 
ity of the barriers the Green functions deviate from the 
above equilibrium values and should be determined from 
Eqs. H]) supplemented by appropriate boundary condi- 
tions describing electron transfer across metallic inter- 
faces. For diffusive superconductors one finds^ 



Aia^GAGi - AiasGsd^Gs 

e^Y- 2Ti,„[Gi,Gs] 

TT V4 + ri,„({Gi,Gs}-2) 
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for the first interface and similarly for the second one. 
Here Ai^2 are the barrier cross sections and <Js,i.2 are 
Drude conductivities of S- and N-terminals. 

Having derived the Green-Keldysh functions G one can 
easily evaluate the current density j in our system with 
the aid of the standard relation 



j - -g^ / tr[f3(GVG)^]d£. 



(6) 



Non-local spectral conductance. The above general for- 
malism enables one to describe electron transport at arbi- 
trary barrier transmissions Ti_„ and T2^„. Here we only 
assume that both NS barriers are sufficiently small to 
provide R12 ^ r = max(rL, tat^ , Tatj). This condition 
allows to effectively linearize Eqs. ([T]) and express the 
solution of linearized Usadel ee^uations via the diffuson 
(lu) and the Cooperon C'' (uj). The diffuson satis- 
fies the following diffusion equation 



DV^]V'- (ij) = S{r~r'), (7) 
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while Cooperon is the solution of Eq. ^ with effec- 
tive charge imbalance relaxation time rg. replaced by 
dephasing time t^. At T ~ tq* depends on the elec- 
tron inelastic relaxation time Tin as^ tq* ~ TmT/ A{T). 

Let us employ the standard representation of the 
Keldysh function G^ = G%-hG'^ with h = fd + hh, 
where /l and /t are respectively symmetric and anti- 
symmetric in energy parts of the distribution function. 
Combining the above expression for G^ with Eq. ([6]) we 
define the current across the first barrier 
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The spectral conductance (71(e) is expressed via the func- 
tions G^, G^. Solving Eq. (P) for G^ and keeping terms 
up to the first order in r/Ri.2, we find 
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where W{e) = in ^ WA^ - for |e| < A, W{e) = 
signe for |e| > A, and 5f™(e) is defined by the stan- 
dard expression^ 
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Note that the terms cx Cs;Ci in Eq. ^ are evaluated 
in the limit Ti^„,T2.n *C 1 where they only matter as 
compared to gf™(e) provided i?i_2 ^ r. The Cooperon 
term oc Ci describes enhancement of Andreev conduc- 
tance by electron interference in diffusive N-metal^i^i^, 
while the term oc Cs accounts for broadening of the den- 
sity of states in the superconductor. The spectral con- 
ductance 52 (s) is given by Eq. ^ with interchanged 
indices 1^2. 

Our next step is to solve the kinetic equation for the 
distribution function fx- In the limit r/Ri2 — » this 

solution is trivial: /•^(e) = and f^^ (e, r^) — h{e, Vj) = 
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(tanh[(£ + eV,)/2r]-tanh[(£- 6^^0/27']) /2 (j = 1,2). 
In the first order in r/i?i 2 the function is determined 
from the diffusion equation 
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where K{e) = d{A - \e\)A^/n^ - e{\e\ - A)e^/n^, and 
fl — 6{A — |e|)il . Resolving Eq. pT|) and substituting 
the result into Eq. ([8]), we obtain 

h{Vi,V2)^ I de[5ii(£)M£,^i)-ffi2(e)/i(£,V^2)],(12) 



where 
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and similarly for the current 12- The last two terms in 
the local conductance 511(e) (|13p describe partial con- 
ductance suppression respectively due to local charge 
imbalance inside the superconductor and due to non- 
equilibrium quasiparticles in the normal metal. At ener- 
gies |e| > A Eq. (fT4|) accounts for the effect of non-local 
charge imbalance which yields non-zero contribution to 
512(e) already in the lowest order in l/i?ii?2. In con- 
trast, at subgap energies this lowest order contribution 
vanishes identically manifesting the well known cancella- 
tion between EC and CAR termsi^. This cancellation is 
lifted in higher orders in barrier transmissions^^. Accord- 
ingly, the full expression for 512(e) (fT4|) docs not vanish 
also for |e| < A and describes non-trivial interplay be- 
tween CAR and direct electron transfer in the presence 
of disorder. 

Eq. (fT4|) for the non-local spectral conductance - to- 
gether with Eqs. ([9]), (fTO|) and (fT3|) - is the central result 
of this work. Note that this result is not specific to par- 
ticular geometry of Fig. 1 but applies for other diffusive 
NSN structures as well. 

Here, however, we will only analyze the system with 
effectively quasi-one-dimensional superconducting and 
normal wires, as shown in Fig. 1. Assuming X2 > xi 
we obtain 
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S,D,^{^+1/t^)/D, 



(15) 



where j — 1,2. Here 5*5,1^2 and Ds,i,2 are respectively 
effective cross sections and diffusion coefficients of the 
corresponding terminals and k — ^J{—iuj -)- \/tq*)/ Ds- 
Substituting Eqs. HH]) into ([T31)-(IT1]) we arrive at the 
conductance matrix describing the system in Fig. 1. 
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FIG. 2: Non-local spectral conductance (712(e) (normalized by 
Go = rxirL-x2/^LRiR2) for diflfusive NSN structures. We set 
L = 10\/25s(0), xi = 0.45L, and X2 = 0.55L. (a) The case of 
two identical barriers with resistances R\ — R2 — -R/e^NchT 
{Nch is the number of channels and r is the barrier transmis- 
sion) and for rjv^ — riVj ~ 0. (b) The case of two tunnel bar- 
riers with Ti,„,T2^„ < 1 and for Lt = L2 = Ljv, r^^ = r^^, 
Ri = R2, tni^n/LiRi = 0.0025 and V2rL^s/LRi = 0.005. 



Zero-bias anomaly. Let us first analyze the tunnel- 
ing limit Ti^„,r2.n ^ 1. In this case at subgap en- 
ergies the term 5f™(e) (fTO|) can be neglected and for 
i?i,2 = Di,2/L\ 2 ^ |e| < A we obtain 
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Here r^^i^e) = rL£,s{£)/L and r^^ ^ie) = rAr^ 25i_2(e)/ii,2 
are Drude resistances of the segme nts of S- and N-metals 
with respective lengths ^s(e) = y/Ds/2i} and ^1,2 (e) = 

At small energies the local spectral conductance di- 
verges as 511(e) (X l/y/e which is just well known 
disorder-induced zero-bias anomaly^i^i^. For the non- 
local conductance P7| this divergence turns out to be 
even stronger, 512(e) cx 1/e, since quantum interference 
in both diffusive normal metals simultaneously enhances 
non-local electron transport in our system. Thus, we 
predict a sharp low energy peak in the non-local con- 
ductance which occurs in the presence of disorder in the 
N-terminals, see also Fig. 2. Accordingly, the differ- 
ential conductance Gi2{V2,T) = —dIi/dV2 increases as 
G12 (X- l/max(eV2,r) with decreasing voltage and tem- 
perature. 

Eq. applies down to e ~ i^i and for even smaller 
energies r^j (e) should be substituted by 2rMi ■ Then for 
TNt , > ?-Cs = (0) we get 



512(0) = Gi2(0,0) 



O r?2 n2 



(18) 



We also note that in the case of strongly asymmetric 
barriers R2 <C -Ri the dominating contribution to 512 
scales as c>c l/i?ii?2 rather than (x l/i?fi?2- 

Turning to the case of high barrier transmissions 
Ti^m 72,Ti ^ 1 we observe that in this case 51 is domi- 
nated by gf^^ (|10p while other contributions can be ne- 
glected. In particular, for fully open barriers at subgap 
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FIG. 3: (a) Non-local resistance Ri2{T) (|20ll . normalized by 
-Ro = r^^TL-x^lTL, at different barrier transmissions r and 
at TQ*,T^ 00. We chose Ri = R2 — -K/e^NchT = h/20e^ 
and set L = 20\/2Cs(0), a;i = 9^Cs(0), = llV2Cs(0), 
VL = 200 r2 and vn^ = rjv^ = 0. (b) Ri2{T)/Ro at different 
|s2 — xi\. Other parameters are the same as in Fig. 2b. 



energies we obtain gi^2 = 9i2 — 2/-Ri,2 and, hence 
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Non-local resistance and charge imbalance peak. Let us 
now define non-local linear resistance 
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Combining this equation with Eq. HH), at T ^ A we 
arrive at a very simple and universal formula 



i?i2 = (r5,/2)e-l^=-^^l/«^ 



(21) 



It is remarkable that independently of both barrier and 
N-terminal parameters the subgap non-local resistance is 
set only by the normal state resistance r^^ of the super- 
conducting wire segment of length and by the dis- 
tance between the barriers measured in units of ^5. At 



low T the dependence i?i^2 ~ roexp[—\x2 — xi\/^s{0)] 
was observed in experiments^ with tq ~ 0.56 Q,. For 
the parametersi^ we estimate ro = r^^ /2 in the range of 
one rj. A similarly good agreement is found between Eq. 
([2T|) and experimental results^. 

The temperature dependence of i?i2(T) is depicted in 
Fig. 3. In the tunneling limit it exhibits a well pro- 
nounced peak which originates from the competition be- 
tween charge imbalance and Andreev reflection2£i2i. The 
maximum value of the non-local resistance i?i2 is reached 
at T* ~ 2A/ln(i?ii?2/r-|g) and reads 
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where A = aL and a — 1 — rj.-^/ri^. Thus, the peak 
resistance i?i2(T*) decreases linearly with increasing dis- 
tance |a;2 — xi| between the barriers. This behavior agrees 
well with recent observations^^. Furthermore, with the 
parameterSiiS we estimate A to be of order a micron in 
agreement with experimental findings. Such values of A 
appear lower than typical values of the charge imbalance 
relaxation length Aq* = \/Dstq' ■ The latter length 
scale is expected to gain importance only for Aq. < A. 

In summary, we developed a microscopic theory of non- 
local electron transport in diffusive NSN systems which 
accounts for non-trivial interplay between crossed An- 
dreev reflection, disorder, quantum interference and non- 
local charge imbalance. Our results can be directly used 
for quantitative analysis of future experiments. 
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